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This talk is about the available hadroproduction data to be used for the tuning of neutrino beam- 
line simulations, going from conventional neutrino beams to superbeams and neutrino factories. 
Comparisons of data with existing Monte Carlo computations will be shown. Recent results from 
the HARP experiment at CERN PS will be summarized, including data useful for a better under- 
standing of extended air showers and atmospheric neutrinos simulations. 
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1. Introduction 

The simulation of a conventional beam with a Monte Carlo (MC) is a delicate task due 
to complicate cascade processes involved in the neutrino production. The paucity of available 
hadroproduction data, needed for MC tuning, can limit systematically the precision in the calcula- 
tions. In addition, new hadroproduction data at low energies are of great interest for extended air 
showers (EAS) and atmospheric neutrinos simulations and for the neutrino factory (NF) design. 
One relevant point is how existing MC simulations compare to available hadroproduction data. 
The more recent dedicated hadroproduction experiments in the field are NA56/SPY [fj, MIPP [|], 
NA61 flf] and HARP [@. At low energies (< 15 GeV), the main experimental results come from 
the HARP experiment at CERN PS and will be briefly summarized here. At higher energies we 
refer to ^ for further details. 

2. The HARP experiment at CERN PS 



The HARP experiment at CERN PS [|J] was designed to study hadroproduction on nuclear 
targets (from H2 to Ta) in the incident momentum range between 3 and 15 GeV/c. The HARP de- 
tector is shown in Fig. [I] and includes different subdetectors for tracking and particle identification 
(PID) over the full solid angle. At large angle (20° < 6 < 160°) tracking and PID are performed by 
a TPC and an array of RPC counters. In the forward direction (d < 14.3° ) the tracking device is a 
set of drift chambers, recovered from the previous NOMAD experiment, while the PID is provided 
by a threshold Cerenkov counter, a time of flight wall (TOFW) and an e.m. calorimeter [Q]. Beam 
particles are tagged by a system of beam TOF detectors (TOFA,TOFB) and Cherenkov counters. 
Data were taken in 2001 and 2002, for a total of about 420 M triggers in ~ 300 experimental 
settings Drift Cnambers 



TPC and RPCs in 
solenoidal magnet! 




Figure 1: Layout of the HARP experiment at CERN PS. 



3. Results for simulation of NF beams. 

The baseline option for a NF target is a Hg jet target with impinging particles at energies 
10 ±5 GeV. Available data are very scarce and for the MC tuning the HARP data on heavy targets, 
such as Ta or Pb, are of utmost importance. The kinematical coverage of the HARP experiment is 
compared with the acceptance of a typical NF design in figure ||. The experiment covers the full 
momentum range of interest for production angles bigger than 0.35 rad. The pion yield increases 
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Figure 2: Left panel: kinematic region in the p — 6 plane covered by HARP as compared to the acceptance 
of the input stage of typical NF designs. Right panel: n + (closed symbols) and n~ (open symbols) yields 
as a function of the incident proton beam momentum for different design of the NF focussing stage. The 
circles indicate the integral over the full HARP acceptance, the squares are integrated over 0.35 rad < 9 < 
0.95 rad, while the diamonds require in addition the momentum cut 250 MeV/c < p < 500 MeV/c. 

linearly with momentum and has an optimum between 5 GeV/c and 8 GeV/c, as can be seen in the 
right panel of figure ||[ Final results for pion production on heavy targets have been published in 
reference [Q] and some comparisons with available MC simulations are outlined in figure [3[ None 
of the considered models describe fully HARP data. However, 7T + production is described better 
than n~ production. At lower (higher) energies binary and Bertini models from GEANT4 (the FTP 
model from GEANT4 and MARS) seem more appropriate. Parametrized models (such as LHEP 
from GEANT4) show relevant discrepancies, up to a factor 3. 

4. Results for simulation of EAS and atmospheric neutrinos 

Results on cryogenic targets, such as N2 and O2 have a direct impact on the precise calculation 
of atmospheric neutrino fluxes and on the improved reliability of extensive air shower simulations 
by reducing the uncertainties of hadronic interaction models in the low energy range. In particular, 
the common hypothesis that p-C data can be used to predict the p-N2 and P-O2 pion produc- 
tion cross-sections may be tested. HARP has published results ||8|] on charged pion production 
cross-sections in interactions of 12 GeV/c protons on C, O2 and N2 thin targets, in the kinematic 
range 0.5 GeV/c < p n < 8 GeV/c and 50 mrad < 6 K < 250 mrad. Some results, showing also a 
comparison with available simulations, are reported in figure 0. 

5. Results for simulation of conventional neutrino beams: the MiniBooNE and K2K 
physics cases. 

Prediction of the far detector spectrum in the absence of oscillations is a key ingredient in a 
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Figure 3: Experimental results from HARP at 12 GeV/c for p-Ta cross sections for 7Z ± production, as 
compared to MC models. See [f7j for further details and more data-MC comparisons. 
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Figure 4: Left: p-C>2,p-A^2,p-C cross sections at 12 GeV/c; right: comparison of n + production in p-A^ 
interactions with different MC models from GEANT4. 



neutrino oscillation experiment. This can be done by an extrapolation from a near detector via a 
nominal far/near ratio estimated by a beamline MC simulation. The error on the observed num- 
ber of events in the K2K far detector (Super Kamiokande) was dominated by contributions from 
uncertainties of normalization (±5%) and far/near ratio (±5%). 

HARP has reported measurements of the 7l + production in p-Al interactions at 12.9 GeV/c 
[Q]. These results have contributed in a significant way to reduce the systematic error associated to 
the FAR/NEAR ratio, thus increasing the K2K sensitivity to oscillation signals [1C]. 
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Similar results were obtained in 8.9 GeV/c p-Be interactions [|Il]] and have contributed to a 
better understanding of the MiniBooNE and SciBooNE v fluxes. Figure || reports the comparison 
with some available MC models. 
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Figure 5: Experimental results from HARP at 8.9 GeV/c for p-Be cross sections for n + production, as 
compared to MC models from GEANT4 (QGSCQGSP) and MARS. 
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